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Abstract 

We propose a generalized extreme shock model with a possibly increas- 
ing failure threshold. 

While standard models assume that the crucial threshold for the system 
may only decrease over time, because of weakening shocks and obsoles- 
cence, we assume that, especially at the beginning of the system's life, 
some strengthening shocks may increase the system tolerance to large 
shock. This is for example the case of turbines' running-in in the field of 
engineering. 

On the basis of parametric assumptions, we provide theoretical results 
and derive some exact and asymptotic univariate and multivariate distri- 
butions for the model. 

In the last part of the paper we show how to link this new model to some 
nonparametric approaches proposed in the literature. 
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1 Introduction 



The setup in extreme shock models is a family {(Xk, Yk), k > 0} of independent 
identically distributed (i.i.d.) two-dimensional random vectors, with Xk the 
magnitude of the fcth shock and Yk the time between the (k — 1) th and the fcth 
shock. The main object of interest is the lifetime/failure time of the system by 
assuming certain schemes for the failure. These models are motivated by the 
possible breakdown of a material or of a system subject to random shocks of 
random magnitude, as it occurs in engineering. Anyway, it is easy to see useful 
applications also in other fields, such as economics, medicine and biology. 
Cumulative shock models and extreme shock models are discussed, as well as 
mixtures of both models, in Gut and Hiisler (1999, 2005) and the references 
therein. 

In the cumulative shock model we consider 

T n = ^2 Y k and S n = 2J X k , 

k<n k<n 

for ri > 1, with Tq = So = 0. The failure of the system occurs if S n > a for some 
n and a. Here a denotes the critical threshold of the system. The time until 
the system fails the first time or the failure time T T with r = min{n : S n > a} 
are then of interest. For results see Gut (1990). 

In the simple extreme shock model one large or extreme shock, larger than a 
given failure (or crucial) threshold 7, may cause the default of the system. The 
lifetime of the system is in this case defined as T v where 

v = min{n : X n > a}. (1) 

This model was dealt with in Gut and Hiisler (1999), and Gut (2001). 
Gut and Hiisler (2005) extended this simple model to a more realistic framework 
by assuming that the failure threshold is not constant, but that it may vary with 
time, depending on the experienced shocks. In detail they assume that large 
but not fatal shock may effect system's tolerance to subsequent shocks, because 
of cracks in the structure for example. To be more exact, for a fixed ao > a 
shock Xi can damage the system if it is larger than a certain boundary value 
P < aJ3. As long as Xi < ao the system does not fail. The crucial hypothesis is 
the following: if a first nonfatal shock comes with values in [/3, ao] the maximum 
load limit of the system is no more ao but decreases to ai S [P, cti]. At this 
point, if another large but not too strong shock occurs in [/3, ao], the new crucial 
threshold is lowered again to «2 G [P, £*i] and so on until the system fails. We 
could call all this "risky threshold mechanism" . Naturally, Vi 

ao > ai > a 2 > ■■■ > P- (2) 

Hence one can define the stopping time v = min \n : X n > otL( n -i)\ with 
L(n) = Y^7=i ^{Xi>p} and L(0) = 0. Gut and Hiisler (2005) have shown that 

1 The value /3 can also vary over time. The only requirement is that it is always lower than 
the corresponding failure threshold. 
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the results for generalized extreme shock models (GEMS) are identical to the 
simple extreme case for nonrandom ctk, while this is not true in the random 
case. We refer to the original paper for more details. 

Even if the modeling of GESM is surely sensible, sometimes it can be worth to 
consider a default threshold a that might even increase, at least initially, say in 
particular for a running-in period of some equipment. This is for example the 
case of turbines' breaking-in in the field of engineering, but other applications 
can be found in electric networks and biology (e.g. Siphonophora in their grow- 
ing process, Dunn et al., 2005). We present in the following section this more 
general model with some theoretical results, which are based on parametric as- 
sumptions. 

In Section 3 we derive some exact and asymptotic univariate and multivariate 
distributions of the parametric model. 

Generalized extreme shock models can even be studied using nonparametric 
techniques. In Cirillo and Husler (2009) a nonparametric urn-based approach 
to extreme shock models is proposed. In Section 4 we briefly show how the same 
approach could be used to model the increasing threshold. 
Section 5 concludes the paper. 

2 Extreme shock models with a possibly increas- 
ing threshold 

In some applications a system has a running-in period during which the critical 
load can increase and the structure is strengthened because at the beginning the 
loads or shocks are large but non-fatal. This may happen in particular until the 
first damage or crack. After such an event the system can only be weakened. 
Such a pattern can be modeled as follows. 

We let the arrival times Tj of the shocks Xi be, as mentioned, a partial sum of 
i.i.d. inter-arrival times Yj with distribution G. The loads Xi, i > 1, arc an i.i.d. 
sequence of r.v.'s with distribution F. A shock or stroke Xi is strengthening 
the material if Xi £ [7, /3)- At the beginning the material supports a maximal 
load a, the critical threshold. After a strengthening stroke, the maximal load 
becomes larger, say a\ = a + b\ with 61 > 0. This boundary increases with 
each strengthening stroke, inducing boundaries <x, = a + bj , j > 1 with bj 7. 
After the first harmful stroke larger than /3, but smaller than the critical level 
at this time point, the load boundary decreases, because of possible cracks or 
some weakening of the material. If it has reached the level (because of k 
strengthening strokes before the first harmful stroke) , the critical level becomes 
now ctk — Ci, and decreases further by the next harmful, nonfatal strokes to 
ctk — C2, Offc — C3, with Cj 7 (> 0). This is shown in Fig. [TJ There might be an 
upper load limit a* for the cek, as well a lower load limit a* for the a* — q. We 
set that a* > /3. It is convenient to set bo = cq = 0. For notational reasons we 
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define the number N-(n) of weakening shocks Xi,i < n, before the n-th shock 
N-(n) = ^2 l i x i € [P,a + b N+{i} - c N _ {l) ), 

with iV_(0) = and the number N+(n) of strengthening strokes Xi, i < n, 
before the n-th shock and before the first weakening or fatal shock 

N + (n) = 5^1(Xi G [7, = 0), 

with iV+(0) = 0. Note that the critical boundary for Xi is on = a + frjv + (i) — 

CjV_(i)- 

In addition, let W be the index of the first harmful shock larger than /?. It could 
indicate even a fatal shock being larger than the critical boundary at this time 
point. Hence, if such a shock occurs: 

W = inf{i : X 4 > /3} < oo. 

If the set is empty, we set W — oo. Hence, the shock Xi has no impact if 
Xi < 7; it induces a strengthening of the material if i < W and X, £ [7, /3); and 
it is fatal, if Xi > on = a + 6jv + (i) — cn-(%)- Note also, that N+(k) = N + (W) 
for all k > W(< 00). 




Figure 1: Realization of a sequence of shocks with strengthening and weakening 
load limits, depending on the values Xi. Here we have v = 32, W = 16, N + (v) = 
2,JV_(» = 2. 

First we analyze in Section 3 the distribution of the number v of shocks until 
the first fatal shock: 

v = min{ j :Xi>a + &jv+(i) - Civ_(i)}- 

where we use bo = and Co = 0. The time until the fatal shock is thus T„. Its 
distribution depends on the distribution of v and G. 
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We consider the asymptotic behaviour of these random variables by letting the 
parameters a = a(t),/3 = /3(t) and 7 = ~f(t) tend to xf < 00 as i — > 00, where 
Xf denotes the upper endpoint of the distribution F. We assume that xp is a 
continuity point of F . For the limit distributions, certain additional restrictions 
will be imposed also on the bk and Cfc, being also dependent on t. Hence v = v(t) 
will tend to 00 in general, depending on the underlying distribution F. 

3 The distribution of v 

The distribution of v can be derived in this more general model as in the basic 
generalized extreme shock model of Gut and Hiislcr (2005). First we derive the 
exact distribution and then analyze the asymptotic distributions which depend 
on the behaviour of the sequences b k and Ck- For the derivation we use the 
notation 

och,l = a + bk — ci 

for any k, I > 0. 

3.1 The exact distribution 

To derive P{v > m} we have to condition on the other random variables. If 
7V_(m + 1) = 0, then W > m and simply P{v > m, iV_(m + 1) = 0} = F m (/3). 
If I > with k < j and m > j + I — 1, we consider the joint distribution 

P{v > m,N+(m) = k,N_(m + 1) = l,W = j} 



j - 1 
k 



x 



F'- l -\l) [F{ 1 )^F{P)] k 
1 



7_ 3 )F m -i- l +\p) JJ[F(/3) -F(a k ,h-i)} (3) 

' h=l 



or for m > j + I and k < j 
P{v = m,N+(m) = k,N-(m) = l,W = j} 

\F j ~ 1 ~ k {~l) [F( 7 ) - F{P)] k 



3-1 
k 



m — j — 1 
I - 1 



F m-j-l w -Q [P(J3) _ F( akih _i)]F(a kll ) . (4) 



h=l 



If k > j or m < j + I — 1, the latter probabilities are 0. 

By summing the appropriate terms, we get the exact univariate and multivariate 
distributions for v, N+(m), N_(l) and W, as well as for N+(v), N-(u). For 
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instance, the joint distribution of 7V+(i/), N-(u) with I > 1 is 

P{N+(v) = k,N_(v) = 1} = ^P{v = m,N+(m) = fc,iV_(m) = 1} 

ill 

= E ( j fc ( m jli 1 )^- fc - 1 ( 7 )[ J F( 7 ) - F{P)] k F m ^-\P) 



X 



\{l-F{a k>h )/F(f3)}F l {f3)F{a k>l ) 



= E r fc V~ fc ~ 1(7) E rr?; M^'-wca 

J>k+1 ^ ' m>j+i ^ ' 

i-1 

x [F( 7 ) - F((3)] k II { 1 " ^K^) W)}^K0 ( 5 ) 

Note that the sums are summing all the probabilities of a negative binomial 
distribution, hence 

P{N + (u) = k,N-(u) = 1} = J2 ( H t k ) Fh ^ 

h>0 
i-1 



fc 

x [F( 7 ) - J] {1 - F(a k , h )/F(p)}F(a k>l ) 

h=0 
2-1 

= F- k -\<Y)[F(>Y) - F{P)] k - F(a Ml )/F(/3)}F(a fc ,) 

= [1 - F(/3)/F( 7 )] fc f[{l ^K, h )/F(^)}[F(a M )/F( 7 )] (6) 



/i=0 

For the case I = we get in the same way 



P{N + {v) = k,N„{v) = 0} 

= Y] = m, 7V + (m) = k,N-(m) = 0, PF = 77;,} 



7n>k 



[l-F(/3)/F( 7 )] fc [^(a M )/^( 7 )] (7) 



Other exact distributions can be derived in the same way by appropriate sum- 
mation. Sometimes, for the derivation of the asymptotic distributions we have 
to approximate these sums to simplify the formulas. 
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We give for later use the following exact distribution. 

P{v > m,N-{y) > 0} 

E P{v> m,N+(m) = k,N-(m + l) = l,W = j} 

fc>0,Z>0j>l 
k>0,l>0,j>l ^ ' 

'™~Af™-1- 1 +\0) f[[F(P) - F(a k ,h-i)}- (8) 

' h=l 



3.2 The asymptotic distribution 

For the asymptotic behaviour let a(t) — > oo,/3(t) —> oo and also j(t) — > oo 
as t — > oo. The asymptotic behaviour depends also on the assumptions of the 
sequences bk and c k , which may depend also on the parameter t. But mostly 
we do not indicate the dependence on the parameter t. From the above finite 
distributions, it is reasonable to use the conditions F((3(t)) — > as t — > oo, 

and 

F { a(t) + b k (t) c l{ t)) = %|) 

Obviously, the a k ,i are monotone by the assumed monotonicity of the sequences 
bk and Cfc, i.e., monotone decreasing in fc with I fixed, and monotone increasing 
in I with k fixed. We consider only the interesting cases with g,a k ,i 6 (0, 1). 
The simplest case occurs if a k ,i = a £ (0, 1) for all k, I. This implies that 

I 

J] [F(f3) - F(a + b k - Cfc_i)] ~ - a)]' 

/i=i 

for each I > 1. But in general we approximate 

i i 
[] [F(/3) - F(a + b k - c h _!)] ~ 

h=l fc=l 

for each / > 1. We use for notational reason also 17/1=1 (1 — a k,h) = 1- 
To simplify the notation, we do not indicate the dependence on t in the follow- 
ing, e.g., we write v instead of v{t). 
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Theorem 1. If fP|] and U0\) hold with g,a k ,i <G [0, 1], then for any k > and 
I > 

i-i 

P{N+(v) = k, N_ (i/) = 1} -> 5 (1 - «y) fe [] (1 - a fc , h )a fc>i 

as i — > oo. 

fYoo/. Use @ for I > 1 and © for I = 0, with the assumptions © and (JTDJ to 
derive immediately the claim, as t — > oo. □ 

For the particular simple case that a k ,h = a for all A: and fo, this limit 
distribution is the product of two geometric distributions 

P{N+{v) = k,N-(v) =1} -» 5(1 -s) fe (l - 

Because of the particular assumptions, the number of strengthening strokes does 
not have an influence on the number of weakening strokes asymptotically, which 
shows the asymptotic independence of N + (i>) and N_(v) in this case. 

Theorem 2. If (0|) and UP]) hold with a k ,i G (0, 1] for each fc, I. 

i) Then for k > = I 

/>oo 

lim P{v > z/F{P),N+{v) = k,N-{v) = 0} = / v k e~ v dv g (1 - g) k a k /k\ . 

J*ia 

ii) For I > 1 and k > 

lim > z/F(B),N+(v) = k,N-(u) = 1} 

/OO rl 
J - y)^ 1 c^{-yu{g- 1 - l)}dy exp{-w} u k+l du 

i-i 

x((l-g)/g) k H(l-a k!h )a k , l /m-l)\) 

h=0 

as t — oo. 
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Proof. Let k > = I and zp = z/F(P). 

P{v > zp,N + (v) = k,N_{v) = 0} = 

= V P{v = m,N + (m) = k, N-(m) = 0} 



111 > Z 



E 1 )^ m " fc " 1 (7)[^(7)--FC8)]*^(«* ) o) 



E ^^ m (7)[^(7)-^(/3)]^K,o) 

m— 

- / « fc exp{- w F( 7 )/F(/3)} d« (F(f3)) 1 [F( 7 ) - F(/3)] fe F(a M ) 



1 f 00 

77 / y k e- y dy(l- g) k ga kfi 



as i — > oo. If A; = 0, then the integral is simply e~ z /° . 
Now let I > 0. By (g]) we have 

> z p ,N+{v) = k,N_{v) = 1} 
= ^ P{v = m, N + (m) = k, N-(m) = 1} = 

m~>zp 

= e e ( j ~ k i )( m ;i~ i y-^)^- i {0) 

m=Zf3 j=k+l \ / \ / 

(-1 

x^ fc ( 7 )(l - g + o(l)) k F l+1 (P) [] (1 a M + o(l))(a fc>/ + o(l)) 

/i=0 

Split the inner sum into three parts, with j < em, em < j < (1 — e)m and 
j > (1 — e)m where we used W = j in the summands. The first and the 
third sum are asymptotically negligible as e — > 0. This can be shown through 
derivations similar to the following ones for the second sum. In fact, for the 
second sum we have 

P{v > Zf) ,N + {v) = k,N-(u) = I, em < W < (1 - e)m} = 



i - 1 

x jj*( 7 )F ,+1 G8)(l - g + o(l)) k [] (1 a*,/. + o(l))(a feji + o(l)) 



m—Zj3 j>em x ' 

(-1 



/i=0 
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oo [(l-e)m] , . . s. 

V V Q-l)!(m-j-l)! 

^ ^ fc!(i-l-fc)!(Z-l)!(m-i-Z)! WJ 



xF( 7 )F i+1 ( / 5)(l- 9 ) k n(l-^K, ! 

oo [(l-e)m] fc 

E E 4^4) ^- fc - 1 (7)^-'(/?) 

m=z@ j>em 



l-l 

x F fc ( 7 )^ +1 (/3)(1 - g) k a k , h )a k>l =: P t 

h=0 



Now we approximate the sums by integrals 

oo 



Pt ~ J (ym) fc (m(l-j/))'- 1 F am ( 7 )F m ( 1 -^(/3)^(m/fc!(/-l)!) 

l-l 

x J F fe ( 7 )F'+ 1 (/3)(l - 

oo „i_ c 

= E / y fe (i -y)^ 1 ^ ym (7)f m(1_!/) (/3)rf2/TO fe+i 

m=z,3 ^ £ 

x F fc (7)^ +1 (/?)(1 - 9? II^ 1 «mKi/(*K* - 1)0 

'\ fc (l-y) , - 1 f*"/^)( 7 )F( 1 -») ,, / fi WC9)d I /« fc+l du 

(-1 

x F fc ( 7 )^- fe (/3)(l - ,g) fc - ak,h)au/W - 1)0 

/i=0 

x <T fc (l - ,9) fc - a klh )a k ,i/(k\(! - 1)!) 

h=0 

The factor FW^^irti-vW-FtfO^) can be approximated for large t by 
cxp(-(l + o(l))yuF(j)/F(f3) - (1 + o(l))(l - y)«) ~ cxp^M.^ 1 - 1) - «]) 
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II 



uniformly for u bounded. Hence by the dominated convergence we get 

/oo pi — e 
J y k (l - y) 1 - 1 cxpi-yuig- 1 - l)}dy cxp{-u} u k+l du 

l-i 

x ((1 - g)/g) k !](! - ak,h)a k ,i/(kl(l - 1)!) 

h=0 

1 

y k (l - y) 1 ^ 1 cxp{-j/u(,g _1 - l)}dy cxp{-«} u k+l du 
i-l 

x ((1 - 9) hit a*,/»Ki/(ft!(i - 1)0 

as e 0. □ 

If we set z = 0, the integrals of both statements can be determined explicitly, 
which implies the result of Theorem 1 in both cases I = and I > 0. 
Other limit distributions are determined by summing the terms P{v = m, N+(m) 
k,N-(v) = l,W = ]}■ Only under additional assumptions these distributions 
can be simplified. Let us deal with such a particular case which generalizes 
Theorem 5.1 in Gut and Hiisler (2005). The formulas and the sums or inte- 
grals can be simplified for instance, if the impact of the strengthening strokes 
is asymptotically negligible, i.e., when a^ h = a hi f° r a U k,h. In this case we 
consider the limit distribution of v and T„. 



Theorem 3. If (0) and il0\) hold with ak.h = a-h & [0, 1] for each k and h > 1. 
Then for zp = z/F{(3) with z > 

oo i l—l 

P{v> Zfi } -> E7r e " z Il( 1 - a ") 

;=o ' h=o 
= 1-H(z) 

and 

P{T V > z p ] -> 1 - H(z/n) 
as t — > oo, where fi — E{Y\) < oo. 

Proof. The assumptions imply that 0/i=o(l — a k,h) = nL=o(l ~ a h) ='■ Q>i- Let 
d a = 1. 

i) For the first statement we start with ((HJ) and apply the assumptions (|9]) and 



11 



(fTU|) to derive the limiting distribution. 



P{v > m,N-(m) > 0} 

m m—j+1 j—1 , . 

= E E Ef fc > H "wi%)-w 

3=1 i=l k=0 ^ 



fc 



X 



' h=l 



E E Ef/l^wra-m 



i - 

m m—j+1 j — 1 

k 

3=1 i=l fc=0 v 

x(^-fjF m -'- l + 1 ^)F'{fi)a l 

Now, using that the sum on fc can be simplified since it is a binomial sum, we 
get for this sum (F( 7 ) + F(-y) - F(/3)) J_1 = (F(/3)) J_1 . 

m m— 3+1 , .s 

= E E ^- 1 G8)(7_i) ii,m ' , " ,+1 ^^) a « 

3=1 i=l ^ ' 

m m — £+1 / .\ 

= E E (?" 5 )^Wo*)3, 



J=l 3=1 



(=1 
i=i f - 



using the normalization m = z/F(/3), which tends to oo. 

ii) The second statement is immediate by applying the weak law of large numbers 
for T v jv — > fi, i.p., as t — > oo. □ 

In other more general cases we have to sum the terms of Theorem 2 to get 
the limit distribution of v and T v . Notice also that the dependence between 
and Yj~ has no influence on the limit distribution of T„. 

4 A link to urn-based shock models 

A nonparametric approach to shock models has been recently proposed in Cirillo 
and Hiisler (2009, 2010). These nonparametric models are based on combina- 
torial processes, and in particular on combinations of Polya-like urn schemes. 
Extreme shock models are modeled in Cirillo and Hiisler (2010) using a special 
version of the reinforced urn process of Muliere et al. (2000), that allows for a 
Bayesian nonparametric treatment of shock models. 
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Generalized extreme shock models are instead modeled in Cirillo and Hiislcr 
(2009) by the means of a particular triangular Polya-like urn. This new model 
is called urn-based generalized extreme shock model (UbGESM). In the next 
lines we aim to show that a similar construction could also be applied, with some 
modifications, to the increasing threshold shock model we have introduced in 
Section [U 

For completeness, let us briefly recall the UbGESM. 

The basic characteristic of the UbGESM construction is to get around the defini- 
tion of the decreasing threshold mechanism of generalized extreme shock models, 
as developed in Gut and Hiisler (2005), by creating three different risk areas for 
the system (no risk or safe, risky and default), by linking every area to a par- 
ticular color and by working with the probability for the process to enter each 
area. If every time the process enters the risky area the probability of failing 
increases, and this can be obtained with a triangular reinforcement matrix, such 
a modeling can be considered a sort of intuitive approach to generalized extreme 
shock models. In some sense, reinforcing the probability for the system to fail 
is like making the risky threshold move down and vice-versa. 
The authors consider an urn containing balls of three different colors: x (safe), 
y (risky), and w (default). The process evolves as follows. 

1 . At time n a ball is sampled from the urn. The probability of sampling a 
particular ball depends on the urn composition after time n — 1; 

2. According to the color of the sampled ball, the process enters (or remains 
in) one of the three states of risk. For example, if the sampled ball is of 
type x, the process is in a safe state, while it fails if the chosen ball is w; 

3. The urn is then reinforced according to its reinforcement matrix RM (bal- 
anced and constant over time) . It means that if the sampled ball is of type 
x, then (it is replaced and) 9 x-balls are added to the urn, if the sampled 
ball is of type y, then d y-balls and A z-balls are added, and if the sampled 
ball is of type z, then 9 z-balls are added. 



RM = 



x y w 

x r 1+0 o o 

y 1+5 A 
w 1 + 



where A = 9 — S 



(11) 



The distribution and the main properties of the urn process can be described 
analytically through the analysis of its generating functiord, see the details in 
Cirillo and Hiislcr (2009). In particular, we can quickly state the following 
theorem. 

2 The generating function of urn histories is a generating function of the form H(z; x, y, w) = 
E^Lo f"( x ' 2/ 1 W )^T> where f n (x,y,w) is a counting function that counts the number of x, 
y, and to-balls in sampling sequences of length n. The gfuh is then that generating function, 
which enumerates all the possible compositions of the urn at time n, given its reinforcement 
matrix and initial composition. 
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Theorem 4 (Cirillo and Hiisler (2009)). Let X n , Y n and W n represent the 
number of x, y and w balls in the urn at time n. Their moments show to 
be hypergeometric functions, that is finite linear combinations of products and 
quotients of Euler Gamma functions. In particular, the moments of order I are 
given by 



E[{X n ) l ] = e l ^f> M n l + 0(n 1 - 1 ), 

( tp-ap \(0 

E \{W n )i] = \ l 6 * n l * + 0(nV-Vi), 
(t v't) 



where to = ao + bo + cq, A = 9 — 8 and i-)^ = r ^"' ) represents the standard 
Pochhammer formula. 



Other results about the limit law of X n , Y n and W n can be found in Cirillo 
and Hiisler (2009), together with results about the asymptotic exchangeability 
of the triangular urn process and its use from a Bayesian point of view. 
The UbGESM shows to be very flexible and it is able to indirectly reproduce all 
the main results of Gut and Hiisler (2005). For example, computing the prob- 
ability that Yio = bo + 1 is like asking which is the probability for the system 
to overcome the risky threshold for the first time in n = 10. In the same way, 
P [W n = co + 1] represents the probability for the model to fail at time n. 



4.1 An urn model for the increasing threshold 



To model the possibly increasing threshold using the urn-based approach, we 
need to introduce a 4-color (x, u, y, w) urn, with initial composition (a , do, b , Cq). 
A possible reinforcement matrbo can be the following: 



RMo 



x 
u 

IJ 

w 



x u y w 
1+9 
9 10 
1+8 A 
1 + 



(12) 



This matrix tells us that the balls of color x, y, w behave as in the standard 
UbGESM, while balls u represent the strengthening shock. In fact, every time 
a it-ball is sampled, the ball is not replaced in the urn and 9 x-balls are added 



3 We could also think of a sacrificial urn (Flajolet et al., 2006), in which ui-balls are removed 
every time a u-ball is sampled. 
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instead. For what concerns the process, it remains in the state it is actually 
visiting. In this way, some strengthening shocks may increase the probability 
of entering or remaining in the safe state, indirectly reproducing the increasing 
threshold mechanism for running-in. If we want all the u-balls to be removed 
sooner or later, in order to avoid further strengthening shocks, we can for ex- 
ample modify the second row of the RM2 matrix, by changing 9 with 9 + 1 and 
1 with 0. 

Theorem 5. Consider an urn process characterized by the reinforcement matrix 
RM2 and with an initial composition (do, <io, °o> Co) of balls. The 5-variables 
generating function of urn histories is: 

H(z;x,u,y,w) = x a °u do y b °w Co (1 - 9x e uz)-^ L (1 - 9w e z)~^ 

• (l -!/*«;-* (l-tl-fou**)*))"*. 

Proof. First notice that states x, u are not directly dependent from states y, w 
through reinforcement and vice versa. Hence matrix RM2 can be seen as a 
combination of a records urn (Flajolet et al., 2006) for x, u and a triangular urn 
for y,w. The result then comes from a direct application of the isomorphism 
theorem of Flajolet et al. (2005). □ 

Given the generating function of urn histories, we can then study the evolu- 
tion of the balls in the urn as in Cirillo and Hiisler (2009). 

Theorem 6. Let X n , U n , Y n and W n represent the number of x, u, y and w 
balls in the urn at time n. Their moments of order I are given by 



( aa+dg \( l ) 

E[(X n ) l ]=9 l \ s n l + 0(n 1 ^), 

(fr 

E[(U n )i] =d 

( to-ao \® , 

E [(W n M = \ l e -L n l l + 0(n«-»t), 
where tg = ao + bo + cq + do, A = 9 — 5. 

Proof. As said matrix RM2 can be seen as a combination of a records urn for 
x, u and a triangular urn for y, w. In particular notice that the number of u-balls 
does not vary over time. 

For what concerns X n , Y n and W n , we only show the proof for y-balls, since the 
methodology is always the same. 
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Set C„ = yf"wn) ■ Taking derivatives of the multivariate generating function, 
(.it J 

one has 



E[(Y n ){\ = E [Y n {Y n - 1) • • • (Y n - I + 1)] = C„ [z ri 



d l H 
dy l 1 



where [z n ] represents the standard notation for the operation of coefficient ex- 
traction (let f(X) = UiZ% b e a generating function, then [z m ] f(z) — 
[* m ]ESo a>iZ i = a m )- 

With some simple manipulations we have that 



E[(Y n ) l ] = C n 



b 



(tri + W)5 

{l-Oz)—^- + 



(J) 



(l-6z) 



(t + (i-l)9)i 



At this point, noting that for 71 < 72, [z n ] (1 - z)-^ 1 = o{[z n ] (1 - z)~ 72 ), 
we discover that only the first term influences the asymptotic behaviour. So, 
thanks to a coefficient extraction with respect to z, we get the desired result. 

□ 

At this point, as shown in Flajolet ct al. (2006) or in Cirillo and Hiisler 
(2009), one can study all the other properties of the urn process. Anyway 
such a study goes beyond the scope of the present section, whose aim is simply 
to build a bridge between generalized shock models with increasing threshold 
and urn-based shock models. It goes without saying that several different urn 
processes can be used to develop the alternative modeling. 



5 Conclusion 

We have proposed an extension of generalized extreme shock model by intro- 
ducing a possibly increasing failure threshold. 

While standard models assume that the crucial threshold for the system may 
only decrease over time, because of weakening shocks and obsolescence, we have 
assumed that, in particular at the beginning of the system's life, some strength- 
ening shocks may increase the system tolerance to large shock, as it happens in 
running-in phases. This is for example the case of turbines' breaking-in in the 
field of engineering. However other fields of applications are easily identifiable. 
On the basis of parametric assumptions, we have provided theoretical results 
and derived some exact and asymptotic univariate and multivariate distribu- 
tions for the model. 

In the last part of the paper we have also shown how to link the new model 
to some recent urn-based nonparametric approaches proposed in the literature 
(Cirillo and Hiisler, 2009). 

Acknowledgements: The present work has been supported by the Swiss Na- 
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